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Helical anisotropy and magnetoimpedance of CoFeSiB wires under torsional stress 

C. Tannous, J. Gieraltowski and R. Valenzuelaf 

Laboratoire de Magnetisme de Bretagne, CNRS FRE-2697, BP: 809 Brest Cedex, 29285 France 

\Instituto de Investigaciones en Materiales, Universidad Nacional Autonoma de Mexico, 

P.O. Box 70-360, Coyoacan, Mexico D.F., 04510, Mexico 



Recent measurements of the magnetoimpedance (at a fixed frequency of 1 MHz) of cobalt-rich 
wires subjected to torsion stress show an asymmetry as a function of torsion angle stemming from 
residual anisotropies induced during wire fabrication. We interpret these measurements with a 
^^' simple model based on the competition between a circumferential magnetic anisotropy and another 

^ I one induced by torsional and residual stress. This allows extraction of the physical parameters of 

the wire and explains the positive and negative torsion cases. The agreement between theoretical 
and experimental results provides a firm support for the model describing the behaviour of the 
OO ' anisotropy field versus static magnetic field for all torsion angles. 



PACS numbers: 75.50.Kj; 72.15.Gd; 75.30.Gw; 75.80.+q 
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INTRODUCTION 



The giant magneto-impedance effect (GMI) in amorphous ribbons, wires and thin films has become a topic of growing 
interest for a wide variety of prospective applications in storage information technology and sensors possessing high 
sensivity and fast response [li, . 

Magneto-impedance effect (MI) consists in the change of impedance introduced by a low-amplitude ac current lac 

i-rt ' flowing through a magnetic conductor under application of a static magnetic field Hdc (usually parallel to the direction 

Pj of the ac current). The origin of this behavior is related to the relative magnetic permeability fir and the direction 

O of the magnetic anisotropy field Hk- When MI measurements are performed on wires, injecting an ac current in the 

1^1 wire and applying a dc magnetic field H^c along the wire axis, one probes the wire rotational permeability wire that 

[ controls the behaviour of the MI. 
Cn ■ Ordinarily, the MI curve is symmetric with respect to the dc magnetic field iJ^^c with a single (at H^c—O) or a 

double peak (when Hfi^. — ±Hjnax)j the anisotropy field of the wire being considered as the value of H„iax- Applying 
a torsion stress to the wire alters its MI behaviour versus field Hdc in a way such that one observes, for instance, a 
r^ ' single peak only for positive stress and a double peak for negative stress. 
\^ Recent work by Betancourt and Valenzuela (BV) 3] tackled MI measurements of CoFeBSi wires subjected to torsion 

^D stress IQJ, at a fixed frequency of 1 MHz, obtaining the aforementioned asymmetry |j,l3 of the MI profile with respect 
"^ to torsion stress. Interpreting the asymmetry in terms of residual anisotropies induced during wire fabrication [3| 

^^ ' leads to an effect that respectively counterbalances or enhances the wire circumferential anisotropy with positive or 
"trt ■ negative torsional stress. 

Cj I The wires were prepared by the in-rotating water-quenching technique and have the nominal composition 

I ■ (Co94Fe6)72.5Sii2.5Bi5. Their typical dimensions are 10 cm in length and 120 /im for diameter. They possess a 

■ O ' characteristic domain structure dictated by the sign of the magnetostriction coefficient A^ which is expected to have 

C , a small value (|Af | ^ 10^^ in these Cobalt rich alloys). When \f is negative, an axially oriented magnetization core 

exists surrounded by circular domains [l|,|3, whereas in the opposite case, radially oriented magnetization states exist 

such as those observed in Fe rich alloys |5J- After performing several measurements such as Barkhausen jump and 

. !^ ' magnetization reversal measurements |3| , BV favor the negative sign for the Af coefficient and presence of helical 

^^ . anisotropy jlflj. 

?H ' In this paper, a simple model for the anisotropy field, based on the competition between an existing circumferential 

anisotropy and a torsion-induced one with presence of residual torsional stress, is proposed. This model allows the 
extraction of the relevant physical parameters of the wire and explains the observed behavior in the cited work [3| as 
a function of the static magnetic field for all torsion angles. 

THEORY 

In order to reveal the residual stress within the wires, positive and negative torsion stress 10r| are applied and 
the effect is analysed with MI response. Following the work of Favieres et al. |3 on cylindrical CoP amorphous 



multilayers electrodeposited on copper wires, we use a stress-induced easy axis making an arbitrary angle 5t [ifl 
with the circumferential direction (see fig.0. We extend this model to include the presence of residual torsion stress 
induced by fabrication. The total magnetoelastic energy in presence of applied and residual torsional stress (ct,.) is: 

E = Kr^s\v?{e) + ^Af (ff + CT^) sin2(^T - e) (1) 

where Kc is the fabrication induced circumferential anisotropy of the wire and A^ the rotational saturation mag- 
netostriction constant jllj (see fig.l). We are in the linear case where the magnetoelastic energy is proportional to 
stress a. Our energy functional does not include demagnetization nor Zeeman terms originating from the dc field 
Hdc applied along the wire axis. Energy due to the ac circular field Hac created by the injected ac current lac is 
neglected. The anisotropy field as a function of stress is determined by identifying it, experimentally, with the value 
of the dc field Hmax where the MI curve displays a peak. Previously Makhnovskiy et al. [2| introduced explicitly 
these terms in the energy functional with a single helical anisotropy []J| term in contrast with our work dealing with 
several competing (intrinsic, residual and induced) anisotropies. The stress a (on the surface of the wire) is related 
to the torsion angle 5 through: 

a = {GaS)/l (2) 

where G is the shear modulus, a is the radius of the wire and I the length (see fig. 2). In general, a has a radial 
dependence since the wire has an inner core with an axial magnetization surrounded by circular domains [l], |2| when 
A? is small and negative. We neglect such dependence and consider, for simplicity, the value of stress on the surface 
[a| of the wire. Introducing the saturation magnetization of the samples Mg (measured with a vibrating sample 
magnetometer [jl as ^qMs = 0.8 Tesla), the energy may be rewritten as: 

E = {MsI2)[Hk sin2(0) + H„ sm^{e - 5t)] (3) 

This shows that we have a competition between a circumferential anisotropy field Hk = 'iKc/Ms and a torsion- 
stress induced one (called henceforth twist field including the residual stress) H^ = 3Af (cr + ar)/Ms ■ Using a 
Stoner-Wohlfarth approach |12l Il3l | , this amounts to rotate the astroid equation from a set of perpendicular fields 
Hx,Hy: 

Hl'^ + ny^ = hT (4) 

to another by an angle Q* resulting in another astroid. This angle is the same that gives the new equilibrium 
orientation of the magnetization and is simply obtained from the condition minimizing the energy, i.e. \dE j dd\^^Q, — 
0. We obtain: 

tan(2r) = ^-g"(^^y) ^ (5) 

^ ' Hk + Ha cos{26t) ^ ' 

Using this expression of 6* ^ the resulting anisotropy field H* is written as: 



H* = ^Hl + HI + 2HkH„ cos{26t) (6) 

The total energy after this transformation is written as: 

E = {Ms/2)[H* sm^{e - 6*) - H* sm^{e*) + H^ sm^iSr)] (7) 

indicating the presence of helical |10| anisotropy (since we have an angle 5t in general different from and 90 
degrees) resulting from the competing circumferential and torsion (appliedH-residual) induced one. In fact we ought 
to find from the MI experiment the angles bounding the value of 5t- 

The MI peak value versus field and torsion angle should correspond to the anisotropy field H* given by the above 
formula eq. Ejas explored in the next section. 



NUMERICAL PROCEDURE 

The MI maximum value field H„iax is to be fitted to the theoretical expression of H* . The latter depends on 
four parameters, Hk the zero-stress circumferential anisotropy field, 6 — 3Xf/Ms the normalised magnetostriction 
constant, the angle St and the residual torsional stress (Jr- These can be reduced to three by normalising all quantities 
by Hk. 

We use a procedure based on a least squares minimization procedure of the curve H* versus a to the set of 
experimental measurements [ci, J/ilj^^ „ where yi = H*{a]iJi). a. represents the set of parameters 6, (5t and cr^. The 
slope of the curve H* versus a that appears to be linear is used in the fitting. This choice is sufficient for the 
description of the method we use, nonetheless our algorithm allows us to select any coherent set of criteria we choose 
to fit the data. 

Hence the couple of minimum equations for the data points and the slope are: 

1 " 

— > [H*{<y.]ai) ~ yi] minimum (8) 

n ^-^ 
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Se| minimum (9) 
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where Se is the overall experimental slope. The fitting method is based on the Broyden algorithm which is the 
generalization to higher dimension of the one-dimensional secant method |l4| and allow us to determine in principle 
two unknowns out of three b, <Tr and 5t- Broyden method is chosen because it can handle underdetermined problems 
(since we have 2 equations and 3 unknowns). 

We consider two cases: 

• Case A: Absence of residual stress (ct^ = 0), then H^: is determined as Hk = H*{a.; 0) and b, 6t are obtained 
with Broyden algorithm. 

• Case P: Presence of residual stress, then Hk is a normalising parameter and b, ar are evaluated with Broyden 
algorithm. 

The parameter Ga/l= 60.82 MPa/radian is obtained directly by averaging over the torsion angles and the torsion- 
stress as done in [Sj. Using a wire radius of a— 60 //m and a length of 1=10 cm, we obtain G = 100 GPa which is a 
little different from the value used in 3 (60 GPa). 

The fitting shown in figure 3 (corresponding to case P as explained in the next section) allows us to predict the 
anisotropy field H*, in the negative stress torsion case, from the experimental results of MI measurements; that is 
interesting enough, since the single peak is dominant (at H^c—O) and the double peak is barely noticeable in the MI 
curve versus Hdc- As an example, we can infer that the anisotropy field H* is zero for a negative value of stress a ~ 
-150 MPa. 

RESULTS AND DISCUSSION 

The fitting procedure in case A, discriminates among the positive and negative A^ values in favor of the negative 
value (implying presence of circular magnetization pattern in the wire) . The fitting quality (equivalent to a x^ test) 
is 0.094 when b= 1.39 (A/m)/MPa (roughly equal to the slope of H* versus a) the angle 5t is almost 180 degrees 
and Af = 4.68 lO"'^. On the other hand, the fitting quality is 1.8 10"^, with b= 1.70 (A/m)/MPa, and 6t ^ 0°. In 
this case Af = -4.55 10~^ agreeing with the value found by BV 3j. 

In spite of the successful fit, the result is not acceptable for at least two reasons: Firstly, with torsional stress the 
angle 5t should be at least 45 ° (or between 45° and 90 ° as in Favieres et al. |8| and Makhnovskiy et al. u\), 
Secondly the zero-stress anisotropy field Hk {Hk = 320 (A/m) '^ is too large for these wires (For Cobalt rich wires 
Hk is on the order of a few tens of A/m. 

Then we move on to discuss the next case P. Starting with Hk = WO A/m, we perform Broyden minimization obtaining: 
b= 2.13 (A/m)/MPa yielding |Af |= 7.13 10""^ while ar= 142.4 Mpa. For Hk^50 A/m we obtain: b= 2.05 (A/m)/MPa 
yielding |Af^|= 6.87 10~^ while ar= 154.1 Mpa. It is remarkable to notice that the fitting curve tapers off for negative 
stress (see fig. 3) while getting better for smaller values of Hk- For Hk='25 A/m we obtain: b= 2.03 (A/m)/MPa 
yielding |Af |= 6.83 10"^ while ar= 156.5 Mpa. 



The value |Af |= 6.83 10~^ is slightly larger than found previously by BV. However, we agree with BV H interpre- 
tation of the MI asymmetry in terms of a fabrication induced residual torsion stress (~ 150 Mpa) that is enhanced or 
counterbalanced (since H* — for a — —150 Mpa) with an applied positive or negative torsional stress. 

We conclude that our results are similar to those obtained with the A case and agree as well with those obtained 
by BV despite the use in both cases of the linear fit: H* = Hk + ha that allows to extract from the slope of the H* 
versus a the magnetostriction coefficient with the formula |l3 : 



A 



-{iioMs/i){dH*/da) 



(10) 



In order to examine the general case with cos(2(5t) ^ ±1, we observe that when the magnitude of the twist field 
Ha is small with respect to the zero-stress anisotropy field Hk one may expand, to fourth order (for instance), the 
anisotropy field H* with respect to a obtaining: 



H* = Hk\i + ^cos{25t) + \ 



Ha 



Hk 



sin2((5T) - - 



Ha_ 

Hk 



cos(2<5T)sin^((5T) + 0(fT^) 



(11) 



One may also derive a similar expansion for A^ as a function of a |l5l | . Using the above relationship eq. 1101 one 
gets: 



'(pQMsb/3) { cos(2(57 






Ha 



Hk 



cos{2dT) sin^ {6T) + 0{a^ 



(12) 



Previously, this kind of expansion has been performed by several workers and particularly by Zhukov et al. [ig 
who made also a fitting of the expansion coefficients as a function of frequency. This might be included easily in our 
approach by letting Ha vary with frequency. 

Generally speaking, the linear regime H* = Hk ± ba should not be observed when the angle 6t is different from 
and 90 ° (± an integer multiple of 90°). Then the relationship H* = ^Hf^ + H^ + 2HKHa cos{2St) should apply 
and when the twist field is small {Ha << Hk), the above expressions eg. 1111 and eg. 1121 apply 

Experiments versus frequency involving wires with off-wire axis torsion stress conditions ought to be done in order 
to test the general non-linear case (where Ha not necessarily being very small with respect to Hk), and precision in 
measuring H* versus a should be increased (in the positive torsion angle case of BV l3|) in order to be able to locate 
accurately the peak of the MI versus field Hdc and test the above formula. This work is in progress. 



[2] 
[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9] 
[10] 



[11] 

[12] 
[13] 
[14] 



K. Mohri, T. Uchiyama, L. P. Shen, C. M. Cai, L. V. Panina, Y. Honkura and M. Yamamoto, IEEE Trans. Magn. 38 3063 

(2002). 

L. V. Panina, K. Mohri, T. Uchiyama, M. Noda and K. Bushida, IEEE Trans. Magn. 31 1249 (1995). 

I. Betancourt and R. Valenzuela, APL 83, 2022 (2003). 

J.M. Blanco, A. Zhukov and J. Gonzalez, J. Appl. Phys. 87 4813 (2000), see also; J.M. Blanco, A. Zhukov and J. Gonzalez, 

J. Phys. D: Appl. Phys. 32 3140 (1999). 

G.H. Ryu, C.S. Yu, C.G. Kim, I.H. Nahm and S. S. Yoon, J. Appl. Phys. 87 4828 (2000). 

V. Raposo, J.M. Gallego and M. Vazquez Journal of Magnetism and Magnetic Materials 242-245 1435 (2002). 

M. J. Sablik and D. C. Jiles, IEEE Trans. Magn. 35 498 (1999). 

C. Favieres, C. Aroca M. C. Sanchez and V. Madurga, J. Appl. Phys. 87 1889 (2000). 

D.P. Makhnovskiy, L.V. Panina and J. Mapps, Phys. Rev. B63 (14), 144424 (2001). 

Helical anisotropy is generally accounted for by introducing an anisotropy axis making an angle i/)(45° < V' < 90°) with 

the wire axis (see for instance ref. |3). Note that torsion stress is also called helical because it is modeled as tensile and 

compressive stresses acting perpendicularly on the wire surface with an angle of 45 ° with the wire axis Q . 

The saturation magnetostriction constants in the tensile case (A^) and in the torsion case (A^) are, in general, different 

for a crystalline material. In amorphous materials A^ ~ A^. 

E. C. Stoner and E. P. Wohlfarth, Philosophical Transactions of the Royal Society of London A 240, 599 (1948). 

L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, Pergamon, Oxford, p. 195 (1975). 

Numerical Recipes in C: The Art of Scientific Computing, W. H. Press, W. T. Vetterling, S. A. Teukolsky and B. P. 

Flannery, Second Edition, page 389, Cambridge University Press (New- York, 1992). 



[15] K. Narita, J. Yamasaki and H. Fukunaga, IEEE Trans. Magn. 16, 435 (1980). 

[16] A. Zhukov, V. Zhukova, J.M. Blanco, A.F. Cobeno, M.Vazquez and J. Gonzalez, Journal of Magnetism and Magnetic 
Materials 258-259 151 (2003). 




FIG. 1: Geometry of the anisotropy and stress-induced easy axis. The circumferential anisotropy axis induced by fabrication 
is along e^ and the torsion stress- induced easy axis S makes an angle 5t with e^. The zero-stress circular anisotropy field Hk 
is along e^. In the helical anisotropy case, the magnetization M lies in the 6^,62 plane making a non zero angle 9 with the 
circumferential basis vector e^, otherwise, both the magnetization M and anisotropy axis lie in the er,e^ plane. 
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FIG. 2: Applying a torsion stress to a wire moves the point B to B' by angle 5. The deformation is the displacement per unit 
length, i.e. e — BB' /I = aS/l where / is the length of the wire. By Hooke's law, stress and deformation are related by the shear 
modulus G, hence a — Ge — G{aS/l) 
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FIG. 3: Experimental results of the anisotropy field H* versus stress as obtained by BV "31 from the field value of Hdc at MI 
peak value and theoretical lines for variable circular anisotropy field Hk at zero-stress. As Hk is decreased from 100 A/m 
(upper line) to 50 A/m and finally 25 A/m (lowest line) we get a weaker tapering off of the line for negative stress a and better 
fitting of the experimental results. The residual torsion stress is ~ 150 Mpa. 



